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We construct an optimal set of single-particle states for few-electron quantum dots (QDs) using 
the method of natural orbitals (NOs). The NOs include also the effects of the Coulomb repulsion 
between electrons. We find that they agree well with the noniteracting orbitals for GaAs QDs of 
realistic parameters, while the Coulomb interactions only rescale the radius of the NOs compared to 
the noninteracting case. We use NOs to show that four-electron QDs are less susceptible to charge 
noise than their two-electron counterparts. 


I. INTRODUCTION 

Quantum dots (QDs) are zero dimensional quantum 
systems with many characteristic properties of atoms®^ 
for example, an energy shell structure was found using 
transport measurements.!®^ We are interested in the po¬ 
tential to store quantum informatiorP using the electron 
spin of QDs.!^ We consider the electron configurations 
of a gate-defined QD with the lowest energies as a real¬ 
ization of a quantum memory. Specifically we discuss a 
two-electron configuration of a QD, where the low-energy 
properties are described by a singlet state and the triplet 
states. The = 0 subspace is two dimensional and de¬ 
fines a qubit,!^ which we call singlet-triplet qubit (STQ) 
in the following. El 

STQs realize universal quantum computation for an 
array of QDs when the transfer of electrons between 
neighboring QDs is permitted.EHUl Many experiments 
have shown that high-fidelity quantum gates can be con¬ 
structed for STQs encoded using gate-defined QDs.EH®! 
In the following, we only consider GaAs QDs, which have 
weak spin-orbit interactions such that the spin and the 
orbital parts of QD wave functions are decoupled. Since 
we only work with the Sz = 0 subspace, we are able to 
neglect the Zeeman interaction and always consider the 
Sz = 0 triplet without further distinguishing it from other 
triplet states. 

The Coulomb interactions between electrons, and con¬ 
sequently correlation effects for QD electrons, are more 
important compared to electrons bound to atoms because 
the sizes of the QD wave functions are larger than atomic 
wave functions.!^ Early descriptions of QDs assumed that 
Coulomb interactions provide a large energy contribution 
to each electron that is added to a QD, but the magni¬ 
tude of this energy c ontrib ution is independent of the 
electron configuration.EEU Very often such a description 
is insufficient because correlations ca n indu ce novel ef¬ 
fects in the energy spectrum of QDs.EEl Even more, 
it has been pointed out that correlations can provide 
all the necessary gate operations for universal quantum 

comput at ion .EEZI 


Weakly correlated quantum systems can be described 
in a mean-field approximation, as for Hartree-Fock cal¬ 
culations (cf., e.g.. Ref. [3]). Mean-field calculations are 
simple, and also th e wav e functions of multielectron QDs 
can be constructed.EEl Hartree-Fock wave functions are, 
however, restricted to a single Slater determinant, and 
a Hartree-Fock description usually fails to correctly de¬ 
scribe strongly correlated quantum systems. Large corre¬ 
lations are expected for weakly confined QDs.^^Even for 
realistic QD parameters, numerical calculations proved 
the failure of Hartree-Fock calculations. El 

Computationally powerful methods have been devel¬ 
oped to describe the correlation effects of few-electron 
QDs. El Very often a full-configuration interaction (full- 
CI) method is used to analyze QD wave functions .El Pull- 
CI calculations diagonalize the QD Hamiltonian using a 
large subspace of possible QD eigenfunctions. The pre¬ 
dicted eigenfunctions of the QDs are exact in the chosen 
basis, but a highly accurate prediction of the ground state 
wave function requires a large set of basis states. We 
show that interaction effects for QDs can be described 
more efficiently using an optimized basis set for Cl cal¬ 
culations. Lowdin described a method to construct ef¬ 
fective orbital wave functions for an interaction quantum 
system, which are called the natural orbitals (NOs). El It 
is well known in quantum chemistry t hat th ese orbitals 
are an optimal basis for Cl calculations.EEHl 

We show that few-electron QDs can be described very 
efficiently using Lowdin’s NOs. Starting from a full-CI 
analysis of QDs, we derive the NOs for realistic GaAs 
QDs. We find that the NOs with the highest occupations 
have a very transparent interpretation because they re¬ 
semble the noninteracting eigenstates of QDs that have 
been rescaled as a consequence of Coulomb interactions. 
The explicit constructions of NOs will also provide an 
analytical description of few-electron QDs. The NO de¬ 
scription is a comprehensive tool to analyze QDs without 
disregarding correlation effects, and we point out that 
NOs can be co nstructed independently from a full-CI 

calculation.EEH 

Finally, the NO description is applied to study the co- 
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Figure 1. Noninteracting orbital energy spectrum of a har¬ 
monic QD for the confining strength hwQ — 3 meV and out- 
of-plane magnetic fields Bz according to Eq. 0- The eigen¬ 
states are the FD states | n, 1) that are described by the orbital 
quantum number n and the magnetic quantum number 1. For 
Bz = 0, the FD states are grouped into atomic shells. We la¬ 
bel the nondegenerate lowest shell as K shell and the doubly 
degenerate first-excited shell as L shell. Out-of-plane mag¬ 
netic fields lift the degeneracies. The QD parameters from the 
highlighted region protect a four-electron STQ from charge 
noise, as described in the main text. 


herence properties of STQs. We focus on charge noise 
because fluctuating electric potentials directly couple to 
the qubit’s wave functions and are one o f the lim iting ob¬ 
stacles to realizing a quantum memoryj^ef. [35] 
predicted a perfect insensitivity to charge noise of a STQ 
that is encoded using a four-electron QD.SSl The noise 
protection was derived from a shell filling model for non¬ 
interacting QD electrons.!^ We find that the four-electron 
configuration remains noise insensitive even in the pres¬ 
ence of Coulomb interactions because the NOs match 
very closely to noninteracting QD eigenstates. When 
Coulomb interactions are included in the description of 
QD electrons, we find nearly an order of magnitude in¬ 
crease in coherence for a four-electron STQ compared to 
a two-electron STQ. 

The organization of the paper is as follows. Sec. 
introduces a model to describe few-electron QDs using 
the noninteracting eigenfunctions of a QD. This model 
predicts a high protection from charge noise for STQs 
encoded using four-electron QDs. Sec. in analyzes in¬ 
teracting QD electrons using NOs, which confirms the 
validity of the noninteracting shell filling model. We also 
characterize the use of NOs in Cl calculations. Sec. lElde- 
scribes the coherence properties of few electron QDs un¬ 
der charge noise. We find that a four-electron STQ is bet¬ 
ter protected from charge noise than a two-electron STQ, 
as predicted from noninteracting QD electrons. Sec. [V| 
summarizes the findings of the paper. 


II. DESCRIPTION OF NONINTERACTING 
FEW-ELECTRON QUANTUM DOTS 


A starting description of few-electron QDs is obtained 
using the noninteracting eigenstates of a QD. We consider 
the single-particle Hamiltonian 


h{r) = 


h? 2 
- 

2m 
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UJn 




{x^ + y^)-^h 


( 1 ) 


of a harmonic confining potential in two dimensions with 
the strength hujo and r = {x^y)^. The orbital contribu¬ 
tions of a magnetic field are introduced through the vec¬ 
tor potential in the symmetric gauge A = ^ {y^ —x, 0)^. 
The out-of-plane magnetic field appears in the QD 
Hamiltonian of Eq. 0 through the cyclotron frequency 
ujc = It gives an additional contribution to the 

harmonic confinement in the second term of Eq. 0- Ad¬ 
ditionally it distinguishes the motion of clockwise and 
counterclockwise circulating electrons through the orbital 
moment Iz = j {xdy — ydx)- 

The eige nvecto rs of Eq. Q are the Fock-Darwin states 
(FD states)NJ^ 
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with Q = ycjQ + (^)^. L^^{x) are the generalized La- 
guerre polynomials. The FD states are described by the 
orbital quantum number n and the magnetic quantum 
number 

As in atomic physics, 'ipni {r) are grouped into or¬ 
bital shells. The eigenenergies of (^) are Eni = 
(2n + |/| + 1) fiQ — Fig. shows the eigenenergies 

of the FD states as a function of the out-of-plane mag¬ 
netic field Bz for a GaAs QD with a typical confining 
strength, hwo = 3 meV [13], and a QD electron with 
the effective mass m = 0.067me, where me is the elec¬ 
tron mass. In the absence of magnetic fields, a highly 
regular level spectrum is observed, where we label the 
orbitally nondegenerate ground state as K shell and the 
doubly-degenerate first-excited states as L shell. The or¬ 
bital components of the out-of-plane magnetic fields lift 
all the degeneracies in the energy spectrum, and there 
is a crossing of energy levels from different energy shells 
at even higher magnetic fields. In the case of very high 
magnetic fields (several tesla), the quantum Hall regime 
is approached with the well-known interpretation of the 
FD states as Landau levels. 

We study the limit of moderate magnetic fields, where 
all the state degeneracies are lifted, but the L shell has 
not been crossed by a state from a different shell yet 
(cf. the highlighted region in Fig. [^. We quickly re¬ 
view the noise protection criterion from Ref. [39]. In the 
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Figure 2. NOs of the lowest energy configurations with = 0 of a two-electron QD and a four-electron QD of Huo = 3 meV 
and = 0.35 T. In each case, the ground state is a singlet, and the excited state is a triplet. The NOs are ordered by 
the energy expectation values of the single-particle Hamiltonian from Eq. Q: E — {'ll; (ri,... rw) |h (t’i)I (ri,... rw)). We 

interpret all the NOs by rescaled FD states | n,/^ . The dominant NOs (drawn with bold lines) resemble the predictions from 

a shell filling model: (a) the two-electron singlet has an occupation of 89% in the lowest NO. Two additional NOs with 5.0% 
occupations are needed to sufficiently describe the singlet. The two-electron triplet has nearly 99% occupation in the lowest 
two NOs. All NOs are summarized in Table(b) The four-electron configurations have nearly 50% occupations in the ground 
state orbital. We interpret this configuration as a frozen core. The singlet has the dominant remaining occupation (43.1%) 
in the next excited NO; the triplet requires two additional NOs (each with 24.4% occupation). In every case, some additional 
NOs have occupations of a few percent or less; all NOs are tabulated in Table [HI 


noninteracting shell filling model of Ref. m, for a two- 
electron QD the electrons are paired in the same orbital 
ground state | n, /) = | 0,0) in the singlet configuration, 
while the Pauli exclusion principle forbids this configura¬ 
tion for the triplet. Here, only one electron occupies the 
orbital ground state | 0, 0), while the other electron is in 
the first excited state | 0 , 1 ) . In the case of four electrons, 
the low-energy properties are again well described with 
a singlet-triplet configuration. In every case, two of the 
four electrons are paired in a “frozen core” in the K shell, 
and the remaining two electrons fill the two orbitals from 
the L shell. In the singlet configuration, both electrons 
are in | 0,1), but in the triplet configuration one electron 
is placed into | 0,1), and the other one in | 0, — 1). 


Charge noise directly couples to the charge densities of 
the QD orbitals. The coherence of a STQ is lost if the 
ground state orbital has a different charge density from 
the excited orbital.l^ In the two-electron configuration, 
the charge density of the singlet state differs from the 
triplet’s charge density because | 0, 0) and | 0,1) have dif¬ 
ferent charge densities. For the four-electron STQ, only 
the charge densities of | 0,1) and | 0, — 1) distinguish the 
singlet and triplet charge densities because in every case 
two electrons are paired in a frozen core in | 0, 0). But 
since |0,1) and |0, — 1) are complex conjugate of each 
other, both states have the same charge densities. A 
four-electron STQ will therefore be protected from all 
the noise sources that couple weakly to the charge con¬ 
figuration of a QD. A more detailed noise analysis, which 
also includes the effects of Coulomb interactions, will be 


given in Sec. IV 


III. DESCRIPTION OF INTERACTING 
FEW-ELECTRON QUANTUM DOTS 

Correlations are important for few-electron QDs, 
which would call into question any predictions of an or¬ 
bital model for few-electron QDs that is obtained from 
the noninteracting QD wave functions. To discuss the va¬ 
lidity of the noninteracting model from Sec.jTTj we analyze 
here the eigenfunctions of the A-particle QDs that we 
obtain from an exact diagonalization of the QD Hamil¬ 
tonian with Coulomb interactions. We use full-CI calcu¬ 
lations to derive the exact eigenfunctions of interacting 
few-electron QDs. The interacting QD problem is solved 
in the basis of 50 noninteracting eigenstates from Eq. ([^ . 
Appx.j^gives a detailed description of the numerical pro¬ 
cedure. The orbital part of the A-particle wave function 
(^1, '^2, • • • 5 “^at) is used to construct Lowdin’s NOs.l^ 
For this, we calculate the first-order density matrix 

g (r', r) =N j dr2 • ... • dvN ( 3 ) 

X (r', r2,..., vn) ^ (r, r2,..., vn) , 

which describes the charge densities and all the single¬ 
particle properties (cf. Appx. for a description of the 
method of reduced density matrices). Note that the prob¬ 
ability density 

p{r) = Jdr2-...-rN {r,r2, ■ ■ ■ ,rN)f ( 4 ) 

is (up to the normalization factor) directly related to the 
first-order density matrix: p{r) = ^^(r,r). The prob¬ 
ability densities always integrate to 1. 
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Table I. Dominant NOs of a two-electron QD in the (a) sin¬ 
glet and (b) triplet configurations. We interpret the NOs 

I NO) as FD states \n,l\ according to Eq. ([^, where the 


tilde represents a rescaling of the Bohr radius as — y by 
^ = /3as. The overlap of | NO) can be maximized from the 
FD-overlap= |(n,/|NO)|^ to overlapi = |^n,/|NO)| with 
the given scaling factor /3. We use 13 = 1.23 to calculate 
overlap 2 for all the NOs, which allows a comparison of the 
NOs with the highest occupations, (a) The two-electron sin¬ 
glet has dominant occupation in 0, 0) , but the states 0,1 


and I 0,-1^ are also needed to accurately describe the sin¬ 
glet configuration, (b) The two-electron triplet has by far the 
dominant occupation in |o,l^, which resembles 

two orbitals filled with one electron each. 



(a) 

two-electron 

singlet 


state 

occupation FD-overlap 


overlapi 

overlap 2 

T76) 

0.8% 

90.1% 

0.96 

90.9% 

68.7% 

oWi) 

5.0% 

98.4% 

0.92 

99.9% 

85.0% 

m) 

5.0% 

98.4% 

0.92 

99.9% 

85.0% 

576) 

89.0% 

95.1% 

1.23 

99.6% 

99.6% 


(b) two-electron 

triplet 


state 

occupation FD-overlap 


overlapi 

overlap 2 

675) 

0.5% 

97.7% 

0.91 

99.9% 

78.1% 

oWi) 

0.5% 

96.6% 

0.88 

99.9% 

80.8% 

oTT) 

49.4% 

99.1% 

1.07 

100.0% 

96.1% 

oTo) 

49.4% 

98.0% 

1.14 

99.9% 

99.3% 


The NOs are the set of normalized, orthogonal func¬ 
tions {(pa (^)} that diagonalize the spectral decomposi¬ 
tion of the first-order density matrix: 

g{r',r) = '^Q„4>l{r')(f>„{r) . (5) 

a 

Qa are called the occupations numbers of the NOs. The 
interpretation as occupation numbers is supported by the 
sum rule 

We continue with an analysis of the two-electron sin¬ 
glet and triplet configurations of a QD using NOs, as 
in Lowdin’s study of the helium atom.^^l In all the fol¬ 
lowing discussions, the NOs are extracted from a full-CI 
calculation, according to Appx. with the parameters 
hwQ = 3 meV and Bz = 0.35 T. These parameters de¬ 
scribe typical QDs that are used in expe riments to realize 
quantum computation with spin qubitFig. (a) 
shows the NOs with the highest occupation numbers for 
a two-electron QD. We see that the intuitive shell filling 
model is by far more valid in the triplet configuration. In 
this case, nearly 99% of the electronic state occupies the 
two lowest NOs. In the singlet configuration, simple as- 


Table II. The dominant NOs of a four-electron QD in the (a) 
singlet and (b) triplet configurations with the same defini¬ 
tions as in Table 12 {3 — 1.42 eases the comparison of the NOs 
with the highest occupations. In each case, nearly half of the 

occupation is in | 0, 0^ , which is interpreted as a core of two 
electrons, (a) In the singlet configuration, by far the domi¬ 
nant remaining occupation is in | 0,, which agrees with the 
picture of two paired electrons in that orbital, (b) The triplet 
configuration has one electrons in | 0,25%), the other 

one in | 0, —1^ 25%). 


state 


(a) four-electron singlet 
occupation FD-overlap /3 overlapi overlap 2 


w) 

0 . 1 % 

91.7% 

0.97 

92.2% 

29.6% 

M) 

0 . 5 % 

88 . 2 % 

1.07 

91.8% 

49.7% 

0 , 3 ) 

1 . 6 % 

99.1% 

0.95 

99.9% 

55.3% 

0 W 2 ) 

0 . 8 % 

99.1% 

0.95 

100.0% 

62.0% 

T76) 

1 . 0 % 

86.1% 

1.10 

90.6% 

66.6% 

072 ) 

3 . 8 % 

99.8% 

1.02 

100.0% 

72.9% 

oWi) 

4 . 5 % 

95.9% 

1.14 

99.6% 

90.0% 

67T) 

43 . 1 % 

88 . 2 % 

1.26 

99.5% 

96.6% 

o76) 

44 . 2 % 

86.4% 

1.42 

99.4% 

99.4% 


(b) four-electron triplet 


state 

occupation FD-overlap 


overlapi 

overlap 2 

0 W 3 ) 

0 . 3 % 

94.0% 

0.88 

99.9% 

43.9% 

1-1) 

0 . 3 % 

89.9% 

1.06 

91.9% 

46.2% 

M) 

0 . 3 % 

89.9% 

1.06 

91.9% 

46.3% 

0 , 3 ) 

0 . 3 % 

94.1% 

0.88 

99.9% 

43.5% 

OW 2 ) 

2 . 1 % 

100% 

1.00 

100% 

70.0% 

T76) 

1 . 0 % 

86.3% 

1.09 

89.4% 

63.9% 

672 ) 

2 . 1 % 

100% 

1.00 

100% 

70.0% 

oWi) 

24.4% 

89.8% 

1.24 

99.5% 

95.6% 

oTT) 

24.4% 

89.8% 

1.24 

99.5% 

95.6% 

o76) 

44.4% 

86.0% 

1.43 

99.4% 

99.4% 


signment of two electrons to the lowest NO is less valid. 
Only 89% of the electrons occupy the lowest NO; two 
additional NOs need to be added to describe the state 
sufficiently. There is a simple explanation for why the 
description with a single NO is less valid in the case of 
the two-electron singlet compared to the triplet. Because 
the spatial wave function of the singlet state is symmet¬ 
ric, the electrons are always closer together in a doubly 
occupied singlet configuration than for any triplet con¬ 
figuration. The Coulomb repulsion is higher for a doubly 
occupied orbital, and a singlet that is restricted to a sin- 
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gle NO has a high energy. The singlet energy can be 
lowered through virtual excitations, which add small oc¬ 
cupations to excited NOs. 

We find that all NOs still closely correspond to FD 


states 


n,/) , where the Bohr radius as = Pclb is mod¬ 


ified compared to the noninteracting case as 


h 

mujQ ’ 


Note that for a subspace of NOs with equal occupations, 
one has the freedom to choose the basis states such that 
they match the FD states most closely. Table |I] sum¬ 
marizes the dominant NOs of the two-electron configu¬ 


rations. The states 


0,0) and 0, l) 


indeed have the 


highest occupation numbers, and their occupations re¬ 
semble the noninteraction shell fillings of FD states. In 
the singlet configuration, the dominant electron occupa¬ 


tion is in 


0,0 


In the triplet configuration, one electron 


IS m 


0, Oy , and the other one in 0, ly . Note that a scal¬ 
ing factor P > 1 was also used in Ref. [33] to account for 
correlation effects of the helium atom. 


The four-electron STQ should be described by NOs in 
the same way as for the two-electron STQ.I^ We analyze 
the lowest two configurations of a four-electron QD of 
hjjQ = 3 meV and 5;^ = 0.35 T, which are again a singlet 
and a triplet. Fig.[^(b) shows the dominant NOs with its 
occupation numbers. In each case, we expect two elec¬ 
trons to be paired in | 0, 0) in a shell filling model. A wave 
function having this form would have 50% occupation in 
the lowest NO. In reality, the ground state’s occupation 
is lower, but we find that the lowest NOs for the singlet 
and the triplet configurations are very similar, and they 


are rescaled FD states 


0,0) (cf. Table 


0 


In the shell filling model, the occupations of the K 
shell describe the orbital properties of this STQ. For the 
singlet, 50% of the charge should be found in | 0,1), while 
for the triplet 25% of the occupation should be placed 
each into | 0,1) and into | 0, —1). The NOs in Fig. i(b) 
and Table confirm this model with large occupations 


in the rescaled FD states 


0,1) and 


0,-1 


Besides 


these dominantly occupied orbitals, we find occupations 
in exited NOs. Although these occupations are low, they 
are still higher than for the two-electron configurations. 

A general implication for a computational analysis is 
that just a few NOs are sufficient to describe the prop¬ 
erties of few-electron QDs. Even though each electron 
configuration has a distinct set of NOs, the NOs with 
the highest occupation numbers have a transparent in¬ 
terpretation. These orbitals represent the FD states with 
a rescaled Bohr radius. Further, the NOs with the low¬ 
est single particle energies are very similar in the singlet 
and triplet configurations of two-electron QDs. In the 
four-electron configuration, we find an identical frozen 
core for the singlet and triplet, and also the next higher 
NOs agree very well in both spin configurations. Based 
on the NOs that were derived in the previous calcula¬ 
tions, we can compare the energies of all the electron 



Figure 3. Probability densities p(r), according to Eq. 
of the lowest energy states of a two-electron QD and a four- 
electron QD with the confining strength hwQ — 3 meV and 
= 0.35 T. In both cases, the ground state is in a singlet 
configuration and the excited state is in a triplet configura¬ 
tion. The dashed lines show radii where the probability den¬ 
sities drop to 1/e of their maximal values. The two-electron 
singlet has a smaller spread of the charge distribution com¬ 
pared to the triplet configuration. In the four-electron con¬ 
figuration, the probability densities are very similar for the 
singlet and the triplet. All properties can be explained in 
shell-filling description using FD states, as explained in the 
text. 


configurations if we restrict the Cl calculation only to 
the dominant NOs. Table El lists the energies of the 
two-electron and four-electron ground states for full-CI 
calculations and for Cl calculations that are restricted 
to the dominant NOs. When we include the NOs with 
the highest occupations, very accurate descriptions of the 
ground state energies are obtained. For a two-electron 
QD, we can reproduce the results of the full-CI calcula¬ 
tion using only four NOs with the given accuracy. The 
four-electron configurations are reliably described when 
nine NOs are included. 


IV. NOISE ANALYSIS OF FEW-ELECTRON 
QUANTUM DOTS 

We start with a noise analysis for the parameters 
fujjQ = 3 meV and Bz = 0.35 T that were analyzed in 
the previous section. Fig. shows the probability densi¬ 
ties p{r) [cf. Eq. @1 of the two-electron and the four- 
electron singlet and triplet wave functions, which are in¬ 
deed the two states of lowest energies in = 0. In the 
case of the two-electron QD, the singlet and the triplet 
probability densities are distinct from each other, which 
is evident in the different spread of the charge distribu- 
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Table III. Comparison of the energies of the singlet configuration and the triplet configuration for the two-electron QD and the 
four-electron QD at hwo = 3 meV and Bz = 0.35 T (cf. Appx. [A|for a description of the numerical calculations), (a) The energy 
of the two-electron singlet, which is the ground state configuration; (b) the energy of the two-electron triplet state, which is 
the first excited state of a doubly occupied QD. The Cl calculation with the two dominant NOs from Table [I] is insufficient to 
accurately reproduce the results from the full-CI calculation. Using the four dominant NOs from Table gives already a very 
satisfactory description of the energy configuration. We obtain a similar finding for the four-electron configuration in the (c) 
singlet and the (d) triplet configurations. The naive description with two (three) NOs from Table |n] in the singlet (triplet) 
configuration is far away from the full-CI results; but adding more NOs from Table [IT] increases the accuracy of the restricted 
Cl calculations. 


(a) two-electron singlet 

(b) two-electron triplet 

(c) four-electron singlet 

(d) four-electron triplet 

Method 

Energy [meV] 

Method 

Energy [meV] 

Method 

Energy [meV] 

Method 

Energy [meV] 

full-CI 

11.1 

full-CI 

12.1 

full-CI 

40.3 

full-CI 

40.5 

2 NOs 

11.8 

2 NOs 

12.6 

2 NOs 

43.3 

3 NOs 

43.3 

4 NOs 

11.2 

4 NOs 

12.1 

6 NOs 

41.3 

5 NOs 

41.9 

5 NOs 

11.2 

6 NOs 

12.1 

9 NOs 

40.6 

9 NOs 

40.8 


tions. The triplet wave function has a larger spread than 
the singlet wave function, which is caused by the Pauli 
exclusion principle, which forces one electron into an or¬ 
bital excited state for the triplet. In the case of the four- 
electron QD, the singlet and the triplet wave functions 
have very similar probability densities. These findings 
agree well with the predictions from the noninteracting 
orbitals of QDs which are successively filled using a shell 
filling method. 

To characterize the degree of similarity between the 
singlet and the triplet configurations, we introduce the 
distance D between the probability densities of the sin¬ 
glet configuration ps (r) and the triplet configuration 
PT (r): 

D j dr\ps{r) - pTir)\. ( 6 ) 

D = 1 for two nonoverlapping probability densities, and 
D = 0 for two identical probability densities. For hujQ = 
3 meV and Bz = 0.35 T, represented by the probability 
densities in Fig. we get D ^ S • 10“^ for the two- 
electron configuration and D ^ S • 10“^ for the four- 
electron configuration. 

We continue with a quantitative noise analysis for 
STQs at various magnetic fields, similar to the study 
in Ref. [39]. A weak single-particle perturbation h' (r) 
renormalizes the energy difference 5Est between the sin¬ 
gle I S) and the triplet | T) states. The leading-order 
contribution to 5Est is given by the direct coupling of 
h' (r) to I S) and | T) 

N 

6Est = (^*)I 1^' (^»)l ^)l • (7) 

i=l 

All operators h' (r) that act only on the spatial part of 
the wave function and are velocity independent will only 
involve the difference in probability densities of the sin¬ 
glet ps (r) and the triplet state pT {r)\ 

SEst = N J dr \E (r) [ps (r) - pr (r)] |. 




Figure 4. Radial dependency of the singlet’s (It’D, solid 
lines] and the triplet’s [pr (It’D? dashed lines] probability den¬ 
sities. The curves are shifted relatively to each other for dif¬ 
ferent out-of-plane magnetic fields Bz. (a) The probability 
densities of the two-electron QDs are always distinct from 
each other, (b) For a four-electron QD, the probability den¬ 
sities are similar at small magnetic fields, but they differ for 
larger Bz. 


( 8 ) 
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Table IV. Probability densities of the singlet ps (r) and 
the triplet pr {r) are compared by the distance measure 
^ ~ I / 1/^-5 (^) “ Pt (^)I cind by their difference of the x 

component of the quadrupole moments = f dr pp (r). 
(a) The probability densities of two-electron STQs differ sig¬ 
nificantly. (b) For four-electron STQs, the probability densi¬ 
ties are very similar at small Bz. 


(a) two-electron configuration 


B. [Tj 

D 

Q¥ 

- Qt [nm^] 

1 

0 

0.081 


-101 

8.0 

0.35 

0.081 


-101 

8.0 

0.7 

0.080 


-96 

7.8 

1.05 

0.056 


-15 

1.3 

1.4 

0.138 


-138 

13.2 

1.75 

0.053 


6 

0.6 


(b) four-electron configuration 

B. [T] 

D 

Q¥ 

- [nm^] 

\Qs -Qt 

0 

0.008 


32 

0.8 

0.35 

0.008 


27 

0.7 

0.7 

0.044 


136 

3.6 

1.05 

0.053 


-89 

3.2 

1.4 

0.019 


-4 

0.1 

1.75 

0.087 


-148 

5.7 


Fig-i shows the radial dependence of the probabil¬ 
ity densities of the two-electron and the four-electron 
configurations of a QD for different out-of-plane mag¬ 
netic fields Bz. We find that the probability densities of 
the two-electron configurations are always differing more 
strongly from each other than the probability densities of 
the four-electron configurations. Especially for small |r|, 
the two-electron singlet and triplet configurations differ 
significantly. This finding can be understood using an 
orbital description for the QD electrons. In the singlet 
configuration, both electrons are placed into the same or¬ 
bital and a high charge density is expected close to the 
origin. In the triplet configuration, the antisymmetry of 
the wave function favors the electrons to be placed fur¬ 
ther apart from each other. 

In four-electron STQs, the singlet and triplet proba¬ 
bility densities are much more similar to each other than 
in two-electron STQs (cf. Fig. [3. Especially in small 
magnetic fields, ps (1^*1) and pTlf^l) matching very 
closely. Larger Bz confine the electrons more strongly, 
and at some point (> 1 T) levels of higher energy shells 
are lower in energy than | 0, —1) (cf. Fig. I^. As a conse¬ 
quence, we expect that ps (1^*1) and pT (1^1) become more 
distinct in increasing Bz^ which is also found when ana¬ 
lyzing the pro babili ty densities in Fig. 

Charge nois^^^^^ dephases STQs while coupling to the 
charge densities of the wave functions. A QD is sur¬ 
rounded by many charge traps, which fluctuate between 
being empty or occupied. A charge trap couples to the 


charge distribution of the QD and dephases a STQ if the 
singlet’s charge density differs from the triplet’s charge 
density. We can quantify this difference using the dis¬ 
tance measure of the probability densities D according to 
Eq. More quantitatively, we can use the quadrupole 
moments of the probability density = f dr ViVjp^ (r) 
of a state | 'i/j). We approximate the shifts of the singlet- 
triplet energy difference SBst that is caused by one 
charge trap according to Eq. We assume that the 
charge trap is some distance away from the QD and use 
a multipole expansion of the QD’s charge distribution.!^ 
The leading contribution to 6Est is caused by the gra¬ 
dient of the charge trap’s electric field dnSj: 


SEst « ^ E 

Qt) 

(9) 

ij=l 



Ne 

~ \dx^i + dye2\ 


(10) 


Because we consider a charge distribution with rotation 
symmetry in the x-y plane and identical probability den¬ 
sities of ps (r) and pr (r) in the z direction, 
and = Qf? are the only nontrivial components of the 
quadrupole tensor. The distance d between the charge 
trap and the QD determines the prefactor in Eq. (10): 

\d^ei+dye2\ = 4 ^ 33 - 

Table |IV| shows that the probability densities ps (r) 
and pT (r^re always differing strongly for a two-electron 
QD, but they can be very similar for a four-electron QD 
at weak magnetic fields. For the data in Table |IV| Bz = 
0.35 T is the optimal value to protect a STQ from charge 
noise because D and | are both very small. We 

explained this finding by the filling of two NOs that have 
the same charge densities in Sec. [TTj Note that the four- 
electron STQ also seems to be equally protected at zero 
magnetic fields and ai B = 0.35 T. STQs at B = 0 will 
still have poor coherence properties because the singlet 
configuration is doubly degenerate. Besides the electron 
configuration that is equ^alent to two valence electrons 

filled into 


0,1 


[cf. Fig.|(b)], - 


our numerical calculation 
finds a state of identical energy for two valence electrons 
in O;^!^ 


The STQ encoding in four-electron QDs was consid¬ 
ered to be optimal be cause the qubit states have identi¬ 
cal electrical moments.^^^Ellln a multipole expansion, the 
coupling of a c harg e tr^ to the electric quadruple mo¬ 
ment dominates .S3 Fig. [^approximates the energy shifts 
6Est through the coupling to the quadrupole moments 
from a single-site charge trap, as derived in Eq. (1^. We 
find that SEst is nearly an order of magnitude larger 
for the four-electron STQ compared to the two-electron 
STQ. While a noninteracting theory predicts 6Est = 0, 
interacting electrons still show a strong noise protection 
for four-electron STQs. 

















Figure 5. Shift SEst of the energy difference between the 
singlet and the triplet states of a two-electron STQ and a 
four-electron STQ which is caused by a charge trap that is 
the distance d away from the QD center, as approximated 
by Eq. (10). SEst is nearly one order of magnitude larger 
for the four-electron STQ than for the two-electron STQ at 
huuo = 3 meV and Bz = 0.35 T. Recall that SEst = 0 for 
four-electron STQs in a noninteracting theory. 


V. DISCUSSION AND CONCLUSION 

We have shown that the shell filling of orbitals describe 
the electron configurations of few-electron QDs. We ex¬ 
plicitly studied GaAs QDs and constructed Lowdin’s NOs 
from the interacting wave functions. Coulomb interac¬ 
tions rescale the radius of the FD states compared to the 
noninteracting orbitals of a QD. Correlations addition¬ 
ally require a few weakly occupied excited NOs. Never¬ 
theless, quantitatively we find the same shell filling model 
for STQs with and without interactions. The immunity 
to charge noise of four-electron STQs, which was indi¬ 
cated from the noninteracting QD orbitals, stays valid 
for interacting QD electrons of realistic parameters. In 
addition, four-electron STQs at weak out-of-plane mag¬ 
netic fields will be protected from all noise sources that 
weakly couple to the charge densities.!^ 

More generally we conclude that a Cl calculation that 
is restricted to only a few NOs sufficiently describes 
the properties of few-electron QDs. Instead of deriving 
the NOs from a full-CI calculation, many self-consistent 
methods to derive NOs in Cl calculations are known in 
quantum chemistry! ^^ * ^^ * An example is the iterative-NO 
method that was proposed by Bender and Davidson.l^ 
Starting from a guess for the NOs, a small set of Slater 
determinants is created that is used for a Cl calculation. 
The NOs can be extracted from this calculation. Again 
a set of Slater determinants is created from these NOs 
that is used for a Cl calculation. Iterating the extraction 
of NOs and small basis Cl calculations ideally converges 
very rapidly. Note that the guess for the NOs does not 
need to resemble the noninteracting wave functions. 

Especially for systems with a large influence of cor¬ 
relations, a full-CI calculation requires a large set of 


basis states to obtain an accurate description of the 
ground states. The iterative-NO method is one of 
many appr oaches to reduce the basis size in numerical 
calculationsj ^^ l ^^ ^ While we focus only on QDs occupied 
with maximally four electrons and typical confining po¬ 
tentials, it has been shown that it is very difficult to treat 
weaker confinements or slightly more electrons in full-CI 
calculations.!^ A self-consistent NO construction would 
help to obtain an accurate numerical description for these 
QDs, and this method can also describe systems that are 
too big for full-CI calculations. 

Our finding is also relevant to analytical descriptions 
of few-electron QDs. Earlier studies very often used a 
minimal-basis Hubbard model to describe the properties 
of few-electron QDs (cf., e.g.. Refs. |39l55l56| k The cho¬ 
sen basis was motivated by the analysis of the nonin¬ 
teracting eigenfunctions of a QD potential, but it has 
been questionable to which degree this few-band Hub¬ 
bard model represents QD wave functions when correla¬ 
tion effects are included. We saw that realistic QD wave 
functions represent the noninteracting orbital level spec¬ 
trum of QDs to a high degree, and an expansion of the 
Hubbard model using only a few additional orbital levels 
can describe all correlation effects. Coulomb interactions 
also rescale the Bohr radius of the FD states compared to 
the noninteracting case. The perspective of highly accu¬ 
rate numerical calculations, with a transparent analytical 
interpretation, should further stimulate the study of NOs 
for QDs in the future. 
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Appendix A: Computational Methods 


We describe the numerical analysis of A-electron QDs 
which is used in the main part of the paper. The N- 
electron configuration of QDs is described by 


N 

n=J2h (n) + E! 5 ’ 

(bi) 


h{r) 

g{r,r') 


(j^ + eAf 

2m 


mujQ 


{x'^ + y^), 


e 1 

Aireoer \r — r'\ 


(Al) 

(A2) 

(A3) 


h{r) is the single-particle Hamiltonian that consists of 
the kinetic energy and the harmonic confining potential 
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of a QD of the magnitude e > 0 is the electron’s 

charge, and m is the effective mass. For the kinetic en¬ 
ergy, we introduce the orbital effects of the magnetic field 
with the magnetic vector potential A. Because the out- 
of-plane confining potential is very strong, only the out- 
of-plane magnetic field component is important, and we 
introduce its orbital effects using the symmetric gauge 
A = ^ (^, —x,0)^. The single-p artic le Hamiltonian in 


Eq. 0 can be obtained from Eq. ( |A2[ ) after some tr ivial 
rewriting. The Coulomb interaction g(r^r') in Eq. ( |A3| ) 
is introduced between all the electron pairs (labeled by 
(bj))- ^0 is the dielectric constant, and is the relative 
permittivity. 

Eor our calculations we use typical parameters of GaAs 
QDs. We use the effective mass m = 0.067me, where rrie 
is the electron mass, and = 12.7 is the relative per¬ 
mittivity. Typical GaAs QDs have = 3 meV. The 


Hamiltonian is discretized using a grid of Gaussian ba¬ 
sis functions,!^ and the full-CI calculation is done using 
a basis of 50 single-particle eigenstates of the Hamil¬ 
tonian from Eq. ( |A2[ ) with the lowest energies. This 
number of basis states was sufficient for the convergence 
of the full-CI calcul ations in earlier calculations of two- 
electron STQs,l^^^^ and we can achieve similar conver¬ 
gence with the four-electron STQs (an energy cutoff sim¬ 
ilar to Ref. m was used here to restrict the size of the 
Cl basis). 


Appendix B: Method of Reduced Density Matrices 

This section describes how reduced density matrices 
describe correlated quantum systems.!^ We define an N- 
electron Hamiltonian 


N 


n=J 2 h (ri) + Y^g {vi, rj) 

ihj) 


i=l 


i Eii9(ri,rj) 


= J dr (r) h (r) ^ {r) J dr J dr' T”*” (r) T”*" (r') g (r, r') T {r') T (r) 


= E 


i- \ na,ma' 

9 2 -^ 


na,ma 
ka" ,la'' 


(Bl) 


with the single-particle Hamiltonian h (r) and the inter¬ 
action Hamiltonian g {r^r'). {i^j) singles out all the pairs 
of the N particles, and Yl'ij sums over all the particles 
except identical ones. We introduce also the notations 
in second quantization. (r) is the field annihilation 
(creation) operator of a particle at position r. We intro¬ 
duced a basis {\na) } of orthogonal and normalized wave 
functions. Cna is the annihilation (creation) operator of 
a particle in the state \na). = {na \h\ma') and 

dka'^ia'" ~ {na^ma' \g\ ka "are the spectral repre¬ 
sentations of the Hamiltonians in the basis {|ncr) }. 

Lowdin introduced reduced density matrices to de¬ 
scribe Welectron wave functions.!^ We define the first- 
order density matrix: 

Q (r, r') =N j dr 2 • ... • dr^ (B2) 

X -0* (r, r 2 ,..., rw) ^ (r', ^ 2 ,..., r^) • 

g (r, r) dr describes the probability that one particle is 
found in the area dr around the point r. The first-order 


density matrix can be conveniently rewritten in second 
quantization: 


g (r, r') = |t’*' (r) T (r') 


= E 0Za'rnAr)^m.'{r'). 

na,ma' 


(B3) 


'Ip') is the spectral decomposition of 
the first-order density matrix in terms of the wave func¬ 
tions (pna {r) = {r \'na ). 

Similarly, one defines the second-order density matrix: 

TT (r, r', r, r') = N {N — 1) J dr^ • ... • dr^ (B4) 

X 'ip* (r,r',r 3 ,... ,rN)'ip {r,r'.r ^,... ,riv). 

TT {r^r'^r^r') drdr' describes the probability that one 
particle is found in the volume dr around r and one 
particle is found in the volume dr' around r'. Also the 
second-order density matrix can be rewritten in second 
quantization: 
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TT (r, r', r, r') = (r) (r') ^ (r') ^ (r) 

= X] ^D'X'"^na (^) ^ma' i^') ^ka" {r) (j^la'" {v') , 


na,ma 
ka" da"' 


with the spectral decomposition tt 


ka"M"' 


^ha^lna'^lcr'"^ko 


2 p) of the density matrix. 


r 


= of Eq. (B1), 


\ i / na,ma' 

and the two-particle Hamiltonian, 


(B5) 


(B6) 


It is very convenient to calculate expectation values 
using the reduced density matrix notation. For example, 
the expectation values of the single-particle Hamiltonian 


/ \ / / • -N \ \ ^ na,ma' na,mi 

{2^9ihj))= 2^ ^ka"M'"9ka",l 

na,ma' 

ka"M"' 


(B7) 


are now easy to calculate. 
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